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ABSTRACT 

We analyse the spatial distribution within host galaxies and chemical properties of 
the progenitors of Long Gamma Ray Bursts as a function of redshift. By using hydro- 
dynamical cosmological simulations which include star formation, Supernova feedback 
and chemical enrichment and based on the hypothesis of the collapsar model with low 
metallicity, we investigate the progenitors in the range < z < 3. Our results suggest 
that the sites of these phenomena tend to be located in the central regions of the hosts 
at high rcdshifts but move outwards for lower ones. We find that scenarios with low 
metallicity cut-offs best fit current observations. For these scenarios Long Gamma Ray 
Bursts tend to be [Fe/H] poor and show a strong a-enhancement evolution towards 
lower values as redshift decreases. The variation of typical burst sites with redshift 
would imply that they might be tracing different part of galaxies at different redshifts. 
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1 INTRODUCTION 

Long gamma-ray bursts (LGRBs, see the reviews by 
Meszaros 2006; Atteia & Vedrenne 2009) are energetic 
radiation events, lasting between 2 and ~1000 seconds, and 
with photon energies in the range of keV-MeV. Our current 
understanding of these sources indicates that the emission 
is produced during the collapse of massive stars, when the 
recently formed black hole accretes the debris of the stellar 
core. During the accretion, highly collimated ultrarelativis- 
tic jets consisting mainly of an expanding plasma of leptons 
and photons (fireball) are launched, which drill the stellar 
envelope. Internal shocks in the fireball accelerate leptons 
and produce the 7-ray radiation through synchrotron 
and inverse Compton processes. External shocks from the 
interaction of the jets with the interstellar medium produce 
later emission at lower energies, from X-rays to radio (after- 
glow). Optical afterglow spectra allowed the measurement 
of LGRB redshifts (Mctzger ct al. 1997), locating these 
sources at cosmological distances (z ~ 0.01 — 8.2), and 
revealing that their energetics is similar to that of Super- 
novae (SNe). Some LGRBs have indeed been observed to be 
associated to hydrogen-deficient, type Ib/c supernovae (e.g. 
Galama et al. 1998; Hjorth et al. 2003; Wooslcy & Bloom 
2006; Starling ct al. 2011). Afterglows allowed also the 
identification of LGRB host galaxies (HGs), which turned 
out to be mostly low-mass, blue and subfuminous galax- 
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ies with active star formation (Lc Floc'h et al. 2003; 
Christensen, Hjorth & Gorosabel 2004; Prochaska et al. 
2004; Savaglio, Glazebrook & Le Borgne 2009). 

Although the general picture is clear enough, its details 
are still a matter of discussion. Among other unanswered 
questions, the exact nature of the LGRB stellar progeni- 
tors is still being debated. Stellar evolution models provide 
a rough picture of the production of a LGRB in a mas- 
sive star. According to the collapsar model (Woosley 1993; 
MacFadyen & Woosley 1999), LGRBs are produced during 
the collapse of single Wolf-Rayet (WR) stars. WR stars have 
massive cores that may collapse into black holes, and are fast 
rotators, a condition needed to support an accretion disc and 
launch the collimated jets. WRs have also large mass-loss 
rates, needed to lose their hydrogen envelope before collaps- 
ing, that would otherwise brake the LGRB jet. This model 
agrees with the observed association between LGRBs and 
hydrogen deficient SNe. However, WRs large mass-loss rates 
imply large angular momentum losses that would brake their 
cores, which would inhibit the production of the LGRB. To 
overcome this problem, Hirschi, Moynet & Maeder (2005) 
proposed low-metallicity WRs (WOs) as progenitors. WOs 
have lower mass- loss rates, diminishing the braking effect, 
but also preventing the loss of the envelope. Another pos- 
sibility was proposed by Yoon, Langer & Norman (2006). 
According to these authors, low-metallicity, rapidly rotat- 
ing massive stars evolve in a chemically homogeneous way, 
hence burning the hydrogen envelope, instead of losing it. 
Low-metallicity progenitor models are consistent with differ- 
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cnt pieces of evidence. First, the works of Mcynet & Maeder 
(2005) and Gcorgy ct al. (2009) show that the collapse of 
high-metallicity stars produces mainly neutron stars, while 
those of low-metallicity stars form black holes. Second, 
LGRB HGs have been found to be low-metallicity systems 
(Fynbo et al. 2003; Le Floc'h et al. 2003; Christcnsen et al. 
2004; Prochaska et al. 2004; Savaglio et al. 2009). Finally, 
the analysis of the statistical properties of the population of 
LGRBs suggests that their cosmic production rate should 
increase with respect to the cosmic star formation rate 
at high redshift, which could be explained as an effect of 
the low metallicity of the progenitors, combined with the 
cosmic metallicity evolution (Daigne, Rossi & Mochkovitch 
2006; Salvaterra & Chincarini 2007). Another possibility for 
WR to lose the envelope without losing too much angular 
momentum is to be part of binary systems as proposed by 
Fryer & Hegucr (2005). 

Understanding the nature of LGRB progenitors is be- 
yond the interest of only stellar evolution, black hole forma- 
tion, and high energy astrophysics. The visibility of LGRBs 
up to very high redshifts (z > 8), allows their use as tools to 
explore star formation and galaxy evolution in the early Uni- 
verse. On the other hand, observations of the environment 
and HGs of LGRBs could reveal important clues about the 
progenitors of these phenomena. Given that star formation 
shifts outward within a galaxy due to the depletion of gas 
in the central regions as the galaxy evolves, that the inter- 
stellar medium of galaxies is not chemically homogeneous, 
and that the chemical enrichment is affected by variations 
of the star formation rate and the production of different 
types of SNe, it is expected that both the LGRB positions 
within a galaxy and the chemical properties of the environ- 
ment in which LGRBs occur depend on redshift and on the 
metallicity of the LGRB progenitors. 

Using high-precision astrometry, 

Bloom, Kulkarni & Djorgovski (2002) and Blinnikov et al. 
(2005) have measured the positions of ~35 LGRBs 
with respect to the centres of their hosts, supporting 
the collapsar model against the (now disproved) neu- 
tron star merger model. The question of the metallicity 
dependence of LGRB progenitors could also be inves- 
tigated comparing these data with model predictions. 
The chemical abundances of LGRB circumburst and 
HG environments were investigated by several authors 
(Prochaska et al. 2007; Niino, Totani & Kobayashi 2009; 
Savaglio, Glazebrook & Le Borgne 2009; Schady et al. 
2010). However, only in a few cases of low-redshift bursts 
a direct measure of the metallicity of the star-forming 
region that produced the LGRB is available. At inter- 
mediate redshift observers usually measure the mean HG 
metallicity, while at high redshift they must resort to 
GRB-DLA techniques, which give the metallicity of galactic 
clouds intercepting the line of sight to the LGRB, but not 
necessarily associated with the burst itself (Prochaska et al. 
2007; Rau A. et al. 2010). 

In this paper, we use cosmological hydrodynamical sim- 
ulations which include star formation and SN feedback 
to investigate the predictions of different progenitor sce- 
narios regarding the positions of LGRBs and the chem- 
ical abundances of their environment. Since galaxy for- 
mation is a highly non-linear process, cosmological nu- 
merical simulations (Katz & Gunn 1991; Navarro & White 



1993; Mosconi et al. 2001; Springel & Hernquist 2003; 
Scannapieco et al. 2005, 2006) are the best tools to 
investigate these LGRB properties. In the past, this 
method has been used by several authors to investigate 
different aspects of the LGRB environment and HGs. 
Courty, Bjornsson & Gudmundsson (2004) have shown that 
requiring HGs to have high star formation efficiency, 
the observed HG luminosity function can be reproduced. 
Nuza et al. (2007) developed a Monte Carlo simulation to 
synthesize LGRB and HG populations in hydrodynami- 
cal simulations of galaxy formation, in the framework of 
the collapsar model. They have found that a bias to low- 
metallicity progenitors (Z < 0.3Z©) is needed to ex- 
plain the observed properties of HGs. Campisi et al. (2009) 
and Chisari, Tissera & Pellizza (2010) used semi-analytical 
models of galaxy formation to study the properties of HG 
populations. Particularly Chisari et al. (2010) developed a 
new approach to model the detectability of the LGRBs. Both 
teams explored models with mass and metallicity cut-offs 
for LGRB progenitors, finding that models with a metal- 
licity cut-off could explain the HG properties, and hence 
supporting previous claims that LGRBs are biased tracers of 
star formation. However in these semi-analytical models, the 
spatial distribution of individual stellar populations within 
HGs cannot be investigated. The chemical abundances of 
LGRB-DLAs were investigated using numerical simulations 
by Pontzen et al. (2010), finding that the clouds producing 
the absorption lie at galactocentric distances of the order of 
1 kpc. 

In this work, we use cosmological numerical simulations 
of galaxy formation to construct synthetic LGRB popula- 
tions, which allow us to investigate the properties of indi- 
vidual stellar populations within HGs. Our simulations are 
similar to those of (Nuza et al. 2007), but with a higher res- 
olution, and include the effects of the energy feedback of 
SNe into the interstellar medium. The metallicities of each 
stellar populations can be estimated, and used to construct 
different metallicity-dependent scenarios for LGRB produc- 
tion within the collapsar model. As stated by Chisari et al. 
(2010), the detectability of LGRBs and their HGs is an im- 
portant aspect that should not be disregarded for a proper 
comparison with the observed samples, hence we included 
it in our population synthesis in the same way as these au- 
thors. 

This work is organized as follows. In Sections 2 and 3 
we describe the cosmological simulations of galaxy forma- 
tion used, and our LGRB population synthesis models, re- 
spectively. We present our results and compare the to avail- 
able observational data in Section 4. Finally, in Section 5 we 
present our conclusions. 



2 SIMULATIONS 

We analyse hydrodynamical cosmological simulations per- 
formed with a version of GADGET-3 which includes star for- 
mation, metal-dependent cooling, chemical enrichment, mul- 
tiphase gas and Supernova feedback (for further details see 
Scannapieco et al. 2005, 2006) The simulated regions repre- 
sent periodic volumes of 10 Mpc side and are consistent 
with a A-CDM universe with the following cosmological pa- 
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rameters: fi A =0.7, fi m =0.3, f2 b =0.04, a 8 =0.9 and H o = 100 
h km s _1 Mpc -1 where h=0.7. 

The feedback model considers Type II and Type la Su- 
pernovae (SNII and SNIa, respectively). The energy per SN 
event released into the interstellar medium is 0.7 x 10 erg. 
The model assumes that stars with masses greater than 
8Mq end their life as SNII with lifetime ~ 10 6 yr. Lifetimes 
for the progenitors of SNIa are randomly selected in the 
range 0.1 — 1 Gyr. The chemical yields for SNII are given 
by Woosley & Weaver (1995) while those of SNIa corre- 
spond to the W7 model of Thielemann ct al. (1993). Initially 
gas particles are assumed to have primordial abundances of 
Xh=0.76 and XHe=0.24. The chemical algorithm follows the 
enrichment by 12 isotopes: 1 H, 2 He, 12 C, 16 0, 24 Mg, 28 Si, 
56 Fe, 14 N, 20 Ne, 32 S, 40 Ca and 62 Zn (Mosconi ct al. 2001). 
We would like to stress that this model has proven to be suc- 
cessfull at regulating the star formation activity and at driv- 
ing powerful mass-loaded galactic winds without the need 
to introduce mass-depend parameters (Scannapieco et al. 
2008). 

We analyse two simulations: S230 and S320, which 
have been also used by De Rossi ct al. (2010) to study the 
Tully-Fisher relation obtaining very good agreement with 
observations. S230 has initially 2 x 230 3 with dark mat- 
ter masses of 5.93 x 10 6 M@ h~ x and initial gas mass of 
9.1 x 10 5 M Q ft -1 . S320 initially has 2x320 3 with dark matter 
of 2.20X10 6 Mq h' 1 and initial gas mass of 3.4 xlO 5 M© h' 1 . 
S320 was only run to z « 2 because of lack of computational 
power. We use this simulation to assess possible numerical 
resolution problems. 

From the general mass distribution, we select virialized 
structures by using the friends-of-friends technique and then 
identify all substructures within the virial radius by applying 
the SUBFIND algorithm (Springel et al. 2001). We select as 
simulated galaxies those substructures sampled with more 
than 3000 particles. 

De Rossi & Tissera (2010) found that the mass- 
metallicity relation (MZR) of galaxies in these simulations 
differs at low redshifts from that reported by Trcmonti et al. 
(2004) so that galaxies have lower mean metallicity than ob- 
served although the shape of the observed MZR is very well 
reproduced. Because of this, we renormalized the simulated 
abundances to make them consistent with observations. 
For that purpose, we adopted the results of Maiolino et al. 
(2008) who proposed a model to describe the evolution 
of observed MZR which matched available observations at 
z=0.07, 0.7 and 2.2. With this adjustement, our simulated 
MZR reproduce observations at different redshifts. For il- 
lustration purposes, we show our analysis at four redshifts: 
z = 0, 1,2 and 3. 



3 SCENARIOS FOR LGRBS 

To construct synthetic LGRB populations from the stellar 
populations described by the simulations we adopt the col- 
lapsar model, in which LGRB progenitors are massive stars 
possibly with low metallicity. We investigate four scenar- 
ios in which progenitors have a mass greater than a certain 
minimum m m i n . For scenario 1 this is the only condition, 
while for the others a maximum metallicity Z c is assumed 
for the progenitors. The values of m m i n and Z c were taken 



Table 1. Properties of the four scenarios proposed for LGRBs. 
The values of minimal masses where taken from Chisari et al. 
(2010) with the initial mass function of Salpctcr (1955). 



Scenario 


z c 


TYl m [ n 


Scl 




~ 80M Q 


Sc2 


0.6 


~ 54M Q 


Sc3 


0.3 


~ 23M Q 


Sc4 


0.1 


~ 6M 



from Chisari et al. (2010), who derived them by fitting the 
LGRB rate observed by BATSE, and are listed in Table 1. 

We estimate the number of massive stars in each simu- 
lated galaxy at each analysed redshift by assuming a Initial 
Mass Function given by Salpcter (1955). We included all 
stars born within r c = 100 Myr. This time interval is larger 
than the mean lifetime of SNII progenitors but it allow us 
to minimize numerical fluctuations and it is small compared 
to SFR variations. The selected progenitors defined the sce- 
nario 1. For scenarios 2, 3, and 4 we impose a requirement on 
the mean metallicity, considering only new born stars with 
Z < Z c = 0.6,0.3,0.1, respectively. 

Following Chisari et al. (2010), for each stellar popula- 
tion represented by a particle p with mass m„ (p, z) at red- 
shift z satisfying the above requirements, we calculated the 
number of LGRBs produced as the number of stars with 
m > m min , 



Jm t,(m)dm 

N{p, z) = m, (p, z) 



Jo.iM ° K(m)dm 



(1) 



where f(m) is the Initial Mass Function with O.IMq and 
IOOA/q its lower and upper mass cut-offs, respectively. 
The intrinsic LGRB rate for a stellar population in any sce- 
nario is then 



r(p, z) 



N(p,z) 



(2) 



We are interested in computing observable properties of the 
stellar populations selected by LGRB observations, such 
as metallicites, a-elements abundances, and distances to 
their HG centre. As discussed by Chisari ct al. (2010), selec- 
tion effects introduced by observations can be modeled by 
weighting the properties of simulated stellar populations by 
their contribution to the total observed LGRB rate at the 
Earth. We applied the method developed by Chisari et al. 
(2010) to estimate the probability that a certain LGRB 
produced at a given z could be observed at Earth. How- 
ever, there might be other biases introduced by observations 
which are difficult to model because of their dependence 
on sensitivity and spectral bands of the detectors and tele- 
scopes. Particularly, the afterglow observations, on which 
the precise positioning of LGRBs is based, are usually made 
in the optical range and could be affected by dust absop- 
tion, biasing the samples towards low metallicity systems. 
As claimed by Fynbo et al. (2009), about 40 per cent of 
LGRBs might be dust obscured. Dust effects have not been 
included in our scenarios hence, caution should be taken 
when comparing our results with observations. We will point 
out posible dust effects when appropriate. 
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At fixed redshift, the weights depend only on the in- 
trinsic LGRB rate of the corresponding stellar population, 
becoming 



1 — I — 1 — I — 1 — I — 

Sc2(Z=0.6Z Siin ) 



Pdet(j),z) 



r{p,z) 



(•3) 



where the sum extends over all the stellar populations p' 
producing LGRBs at a given redshift. For any observable 
property X(p, z) of these stellar populations, its mean ob- 
served value at z must then be 



{X{z))=Y J P^{jP,z)X(p',z). (4) 



4 RESULTS AND ANALYSIS 

4.1 Spatial distribution of LGRBs 

We first investigate the spatial distribution of LGRBs in 
our scenarios. For this purpose, we calculate the distance 
between the LGRB and the centre of mass of its galaxy b. 
To eliminate the effects produced by the growth of galaxies 
as the structure in the Universe assembles, we normalize b 
by taking the ratio fc/r op t, where r opt is the optical radius of 
the galaxy, defined as the radius encompassing 83 per cent 
of its baryonic mass (Tissera & Dommguez-Tenreiro 1998). 

In the Fig. 1 we present the distribution of b/r opt for 
the LGRBs at different redshifts (z = 0,1,2,3), weighted 
by their detectability as explained in section 3. We observe 
that LGRB progenitors tend to reside in the inner regions 
galaxies at high redshift, and to be progressively located at 
larger distances from the centre as redshift decreases. This is 
consistent with the fact that the main sites of star formation 
shift outwards as time evolves and the galactic structure gets 
assembled in a hierarchical fashion. This effect is stronger in 
our scenarios with higher Z c , because low metallicity pop- 
ulations tend to be formed in the outer regions of galaxies 
which are less enriched since all simulated systems exhibit 
metallicity gradients. 

In Fig. 2, we plot the median values of b/r op t in our 
scenarios as a function of redshift, together with the avail- 
able observations of the LGRBs positions within their hosts 
(Bloom et al. 2002; Blinnikov et al. 2005). These authors 
measure the distance of LGRB to the centre of their hosts, 
projected onto the plane of the sky, and normalised by the 
galaxy half-light radius r^. As these authors point out, this 
normalisation is a crude way of deprojecting the values of 
6. To transform them into b/r opt , we assume that LGRB 
hosts can be modeled by an exponential disc, for which 
th = 0.52r opt . Fig. 2 shows that our scenarios are con- 
sistent with observations, except at very low redshifts in 
which the observed median value of 6/r op t drops abruptly, 
while our scenarios remain almost constant. By analysing 
the LGRBs contributing to the lowest-z point in Fig. 2, we 
find that almost half of them (3 out of 7) have 6/r opt val- 
ues consistent with zero. Interestingly, the hosts of these 
LGRBs show evidence of interaction or close companions. 
Hence, the presence of nuclear star formation activity could 
be explained as triggered by galaxy interactions as suggested 
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Figure 1. Distribution of LGRB distances to the galaxy centres 
normalised by the optical radius in our scenarios, at z = 3 (solid 
lines), 2 (dashed lines), 1 (dotted line), and (dot-dashed lines). 



by observations (Lambas et al. 2003) and numerical simula- 
tions (Barnes & Hernquist 1996; Mihos & Hernquist 1996; 
Tissera 2000; Perez et al. 2006). Then, the discrepancy can 
be attributed to the fact that our simulated galaxy sample 
does not reflect the effects of this mode of star formation at 
low redshift since our analysed galaxies are dominated by 
systems with low gas reservoir (De Rossi & Tissera 2010; 
De Rossi, Tissera & Pedrosa 2011). A further piece of evi- 
dence for this explanation is provided by a recalculation of 
the lowest-z point, excluding the three quoted LGRBs (filled 
circle in Fig. 2). The new point lies within 3<r of our scenar- 
ios, showing a better agreement than the original one. 

The large error bars of the observations, which origi- 
nate in the low number of LGRBs with precise positions, 
prevent us from using a goodness-of-fit estimator to deter- 
mine the scenario that better fits the observations. However 
the fact that the observed values are always higher than 
the predictions of scenarios Scl and Sc2 implies that it is 
very improbable that these scenarios could explain the ob- 
servations. Hence our results suggest that LGRB progenitors 
would have low metallicities (Z < 0.3Zq). In the case that 
dust effects introduce important biases in the impact param- 
eter distribution, the preference for low metallicity progen- 
itors obtained from Fig. 2 would have to be re-considered. 
Observations providing new insights on the location of dark 
GRBs may help to resolve this issue. 

4.2 Chemical abundances 

To analyse the chemical abundances of the LGRB progeni- 
tors in our scenarios, we use the ratio [Fe/H] as a measure 
of the iron abundance, and [Si/Fe] as a measure of the rel- 
ative abundance of a elements to iron. In Fig. 3 we present 
the distribution of [Fe/H] for the LGRB progenitors in our 
scenarios at different redshifts, wheighed by the detectabil- 
ity in the same way as in the previous section. We observe 
that the abundance of iron increases as redshift decreases in 
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Figure 2. Median LGRB distances to the galaxy centres nor- 
malised by the optical radius as a function of redshift, in our 
scenarios (Scl squares, Sc2 circles, Sc3 triangles, Sc4 diamonds) 
compared with median calculated by the observations (stars) of 
Bloom et al. (2002) and Blinnikov et al. (2005). 
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Figure 4. Cumulative histograms for [Si/Fe] for 2 = 3 (solid 
line), 2 (dashed line), 1 (dotted line), and (dot-dashed line) in 
all scenarios proposed. 
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Figure 3. Distribution of the relation beween Fe/H for our four 
scenarios at different redshifts z = 3 (solid line), 2 (dashed line), 
f (dotted line), and (dot-dashed line). 



high redshift only the contribution of the latter to the in- 
terstelar medium enrichment is significant, rendering stel- 
lar populations rich in a elements. As redshift decreases, 
the contribution of SNIa becomes important, decreasing the 
abundance of a elements relative to iron. 

The trend of [Si/Fe] to decrease towards lower redshift 
is also observed in Fig.5 (right panel), where we plot the me- 
dian value of the ratio [Si/Fe] as a function of redshift. In this 
figure we also observe that the values of [Si/Fe] of the LGRB 
progenitors is lower and evolve more strongly with redshift 
in the scenarios where the metallicity cut-off is more restric- 
tive. In the left panel of Fig. 5 we show the median value of 
the ratio [Fe/H] as a function of redshift. We find that the 
median value of [Fe/H] decreases with Z c as expected for 
cut-offs progressively more restrictive in metallicity. These 
results indicate that the metallicity cut-off tends to elimi- 
nate old stellar populations highly enriched by SNII, located 
mainly in the central regions of galaxies and originated in 
first outbreaks of star formation. This interpretation agrees 
with the shift of the normalized impact parameters of LGRB 
progenitors observed in Fig. 1. 



all scenarios. This can be understood in terms of the chemi- 
cal evolution of the interstellar medium. As time evolves SNe 
contribute to the enrichment of the medium with iron, hence 
stellar populations born at low redshifts exhibit higher iron 
abundances. This enrichment is stronger in scenarios Scl 
and Sc2, where the metallicity cut-offs are not so restrictive. 

In Fig. 4 we present the distribution of [Si/Fe] for LGRB 
progenitors. We observe that, for all our scenarios, they ex- 
ibit a higher [Si/Fe] as redshift increases, indicating an en- 
hacement in a elements at high redshifts. This is consis- 
tent with the fact that SNIa and SNII enrich the interstellar 
medium. Due to the fact that SNIa progenitors have life- 
times ~ lGyr while those of SNII live only ~ lOMyr, at 



5 CONCLUSIONS 

Aiming at understanding the relation between LGRBs and 
star formation, we analysed the spatial distribution and 
chemical abundances of stellar populations producing these 
phenomena. We investigated four different scenarios for the 
progenitors of LGRBs based on the collapsar model with 
different metallicity cut-offs. 

We compared the spatial distribution of LGRBs within 
their HGs in our scenarios and with the available observa- 
tions. We found that in all our scenarios LGRB progenitors 
reside on average in the outer regions of their galaxies at low 
redshifts, shifting toward the centre as redshift increases. 
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Figure 5. Evolution of median values of [Fe/H] (right) and [Si/Fc] (left) of LGRB progenitors in our different scenarios. (Scl squares, 
Sc2 circles, Sc3 triangles, Sc4 diamonds). 



Scenarios favouring low metallicity progenitors tend to pro- 
duce LGRBs further out from the central regions than those 
allowing high metallicity progenitors. The confrontation of 
our models with available observations supports scenarios 
with low metallicity cut-offs, in agreement with previous re- 
sults (Nuza et al. 2007; Campisi et al. 2009; Chisari et al. 
2010). Particulary we best reproduce current available ob- 
servations for a model where LGRB progenitors are massive 
stars with Z < 0.3. Further precise LGRB position measure- 
ments would help to confirm these trends. 

Regarding [Fe/H] abundances of the stellar populations 
producing LGRBs, we found that in all our scenarios [Fe/H] 
increases as redshift decreases. This effect is less conspic- 
uous in the scenarios with low metallicity progenitors, as 
in these cases the metallicity cut-off restricts the chemi- 
cal abundances of the stellar populations producing LGRBs. 
The a-enhancement decreases with redshift in all our sce- 
narios, as a result of the different contributions of SNII 
and SNIa. Contrary to the detected trend in [Fe/H], the 
a-enhancement shows a stronger evolution with redshift as 
Z c decreases. As previously discussed, these chemical trends 
can be understood within the context of chemical evolution 
in hierarchical clustering scenarios. 

Considering that the results on the spatial distribu- 
tion of LGRB progenitors favours low-metallicity progeni- 
tor models, one would expect that the iron abundance of 
the stellar populations producing LGRBs remains low at all 
redshifts with little variations ([Fe/H] ~ —1)- On the other 
hand, one would expect that the a-enhancement strongly de- 
creases with redshift (by 0.2 dex between z — 3 and z — 0). 
This means that, if LGRBs are produced by low metallicity 
massive stars, their location will be shifted on average from 
the central regions to the outskirts of galaxies. If LGRBs 
can trace the chemical properties of the interestelar medium, 
they may map different regions of galaxies at different red- 
shifts. A test of these prediction could be set up as further 
dust-corrected measurements of the chemical abundances of 
the stellar populations producing LGRBs become available. 
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